Back vowels
d <- reps %$% reps[phoneme %in% c("UW","OW","AA","IY") &
!(np1_phone %in% c("L","N","M","R","NG")) &
stress %in% c("S1", "S2"),]
p <- ggplot(d, aes(x=dob, y=X0, colour=phoneme)) +
geom_smooth(lwd=3) +
theme(text=element_text(size=20)) +
xlab("Date of birth") +
facet_wrap(~ rep, scales="free_y", shrink=T, ncol=2) +
scale_colour_manual(values=color_palette)
print(p)

d <- reps_spkr %$% reps_spkr[phoneme %in% c("UW","OW","AA","IY") &
stress %in% c("S1", "S2"),]
p <- ggplot(d, aes(x=dob, y=X0, colour=phoneme)) +
geom_smooth(lwd=3) +
theme(text=element_text(size=20)) +
xlab("Date of birth") +
facet_wrap(~ rep, scales="free_y", shrink=T, ncol=2) +
scale_colour_manual(values=color_palette)
print(p)

Front lax vowels
d <- reps %$% reps[phoneme %in% c("IY","IH","EH","AE", "AA") &
!(np1_phone %in% c("L","N","M","R","NG")) &
stress %in% c("S1", "S2"),]
p <- ggplot(d, aes(x=dob, y=X0, colour=phoneme)) +
geom_smooth(lwd=3) +
theme(text=element_text(size=20)) +
xlab("Date of birth") +
facet_wrap(~ rep, scales="free_y", shrink=T, ncol=2) +
scale_colour_manual(values=color_palette)
print(p)

d <- reps_spkr %$% reps_spkr[phoneme %in% c("IY","IH","EH","AE", "AA") &
!(np1_phone %in% c("L","N","M","R","NG")) &
stress %in% c("S1", "S2"),]
p <- ggplot(d, aes(x=dob, y=X0, colour=phoneme)) +
geom_smooth(lwd=3) +
theme(text=element_text(size=20)) +
xlab("Date of birth") +
facet_wrap(~ rep, scales="free_y", shrink=T, ncol=2) +
scale_colour_manual(values=color_palette)
print(p)

d <- reps %$% reps[phoneme %in% c("IY","IH","EH","AE", "AA") &
stress %in% c("S1", "S2"),]
d$np1_son <- d %$% np1_phone %in% c("L","N","M","R","NG")
p <- ggplot(d, aes(x=dob, y=X0, colour=phoneme, linetype=np1_son)) +
geom_smooth(lwd=3) +
theme(text=element_text(size=20)) +
xlab("Date of birth") +
facet_wrap(~ rep, scales="free_y", shrink=T, ncol=2)
print(p)
